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QCD string model is employed to calculate the masses and spin splittings of lowest 

00 ' 

f-*) ' charmonium hybrid states with a magnetic gluon. Relative decay rates into various 

O ■ 

^ ■ 5- and P-wave D-meson pairs are calculated for these hybrids. 

^ ' PACS numbers: 12.39.-x, 13.25.Jx, 14.40. Gx 

00 ■ 

(N : 

, ; I. INTRODUCTION 

p ^. There exist strong arguments in favour of hybrid assignment for the recently observed 

I y(4260) state 1]. Indeed, this F-meson is definitely a vector one, as it is seen in the 

^ ■ initial state radiation process e~^e~ 'yir'^n' J/ip, but its e~^e~ width is too small for a 

^ ! conventional cc vector, and there is no visible decay into DD pairs, in spite of the large 

2 ' phase space available. It is the latter feature that has prompted the hybrid interpretation 

^ . of the y(4260)j^, as the selection rule is established — see, for example, j^, [4, O, y, [t] — 

O ' which forbids the decay of the vector hybrid into a final state. 
00 . 

O ■ Competing models for the y (4260) exist. One is the [cs] — [cs] diquark-antidiquark model 

• ^ ! 8j. On the other hand, the y(4260) is not far from the DDi threshold, where Di is a P 

?H ' wave 1"*" charmed meson, so the y(4260) state could be associated with the opening of a 



new S-wave DDi threshold [9|. In this regard it is important to assess the consequences of 
the hybrid assignment for the Y . 

Hybrids can be considered as bound states of a quark-antiquark pair and a gluon with 
quantum numbers 

p = C = (1) 

for the magnetic gluon {Ig = j), and 

for the electric gluon {Ig = j ± 1), where Ig is the relative angular momentum between the 
qq pair and the gluon, j is the total angular momentum of the gluon, Iqq is the orbital 



2 



momentum in the quark-antiquark subsystem, and Sgq is the spin of the quark-antiquark 
pair. For a magnetic gluon, the lowest states correspond to Igg = 0, with the 1 hybrid 
being a spin-singlet state with respect to the quark spin, while the J~^, J = 0, 1, 2, hybrids 
being spin triplets. These four states are expected to be degenerate in the heavy-quark 
limit, with the degeneracy removed by spin-dependent quark-gluon interactions. In other 
words, if y(4260) is indeed a ccg vector hybrid, three J~~^ hybrid charmonia, including 
the exotic 1 one, should reside not very far from it. Decays of all these states obey the 
above-mentioned selection rule. 

It is the latter feature, which makes these states potentially interesting. Indeed, the 
same quantum numbers can be also achieved with the electric gluon, 1 hybrid being 
a spin-triplet state with respect to the quark spin, and J^^ hybrids being spin-singlets. 
However, hybrids with the electric gluon couple too strongly to two S-wave final-state 
mesons {D^*^D^*^) and, as estimated in Ref. j^l, do not exist as resonances. On the contrary, 
for hybrids with the magnetic gluon, the modes are forbidden, and the lowest 

possible open-charm modes are the ones with an S-wave D^*^-meson and a P-wave Dj- 
meson, with the thresholds being close to the masses of such hybrids. Due to a limited 
phase space, one expects a considerable suppression of the corresponding decay width, so 
that hybrids with the magnetic gluon could manifest themselves as resonant states, and, as 
such, could be of immediate relevance to the charmonium spectroscopy issues. 

In the present paper we calculate spin splittings in the ccg system with a magnetic gluon 
in the framework of the QCD string model based on the Field Correlator Method (FCM) for 
QCD (for a review of the FCM see Ref. P^). In this method, confining dynamics is encoded 
in gluonic field correlators which are responsible for area law asymptotic for the Wilson loop. 
Starting from the Feynman-Schwinger representation for the quark and gluon propagators 
in the external field, one can extract hadronic Green's functions and calculate the spectra. 
The QCD string model corresponds to the limit of a small (vanishing) gluonic correlation 
length — the so-called string limit of QCD. Then the effective string-type Lagrangian of 
a colourless object (meson, baryon, hybrid, and so on) can be derived. The QCD string 



model was successful 
for example, Refs. 



ly applied to calculate spectra and other properties of qq mesons, see. 
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12| . To account for hybrid excitations one populates the QCD 
string with constituent perturbative gluons |l^. This approach was used before in order to 
investigate the properties of hybrids iJ], and the form of the static interquark potentials in 



3 




Q 

FIG. 1: Wilson loop configuration corresponding to the propagation of the hybrid state, 
hybrids was studied in detail and compared with the lattice simulations in Ref. jisl . 

II. THE MODEL 

In the framework of the FCM, hybrid is viewed as a gluon with two fundamental strings 
attached, with the quark and the antiquark at the ends (see Fig. [1]). Thus the starting point 
of our analysis is the effective Lagrangian of such a system: 



2flg 2/ig 2 

•1 , rl 



- ^1^9-^9!^ dpi^l- If - a\fg - fg\ dp2\/^^^-^, (3) 



h = !f| X + (1 - I3i)^g), 

\rq-rg\ 

h = ^4—^, X + (1 - /Js)^^, 
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where r^, fg, and are the quark, the antiquark, and the gluon coordinates, irig and mq 
are the quark and the antiquark current masses (the gluon is massless), dots stand for 
time derivatives, and the minimal surface is approximated by the straight-line string. The 
Lagrangian ([3]) is written in the einbein field form the einbein fields (or simply einbeins) 
/ig, /ig, and iig are introduced to deal with relativistic kinematics. No time derivatives of the 
einbeins enter the Lagrangian, and the corresponding equations of motion yield second-class 
constraints (see Ref. for the details of the constrained systems formalism and for the 
corresponding terminology): 

dL dL dL ^ 

dflq djJig 

In principle, one can eliminate einbeins using Eq. (jl]). However, only the einbein form of the 
Lagrangian provides a meaningful dynamics for a massless gluon. Application of the einbein 



field formalism to the QCD string was suggested in Ref. [18|, and further developments can 



be found in Ref. 



Starting from the Lagrangian ([3]) one can arrive at the Hamiltonian of the qqg system. 



The general procedure outlined in Refs. UJ, llSl] is rather complicated due to the presence of 



square roots in ([3]). However, if one is interested in the low-lying part of the spectrum, the 
potential-type regime can be considered. To this end notice that the angular velocities /f 2 in 
the Lagrangian ([3]) describe the contribution of the proper inertia of the rotating string. For 
the lowest approximation it is sufficient to retain only the first terms (linear confinement) 
in the expansion of the string terms in powers of 3- Corrections to this potential regime 
coming from the further expansion of the string terms are a footprint of the underlying 
string dynamics. The leading correction of this type, of order 35 is known as the string 



correction 



18j — it will be taken into account later. 



Then the zero-order Hamiltonian takes the form: 



-"0 — ^ \ ^ \ ^ H h (T\rg — rg\+ a\rq — Vg] + Kcoui, l^j 

where the long-range confining force is augmented by the short-range Coulomb potential, 

2\rg-rg\ 2\rg-rg\ 6\rg-rg\ 



The coefficients in Eq. (EI) correspond to the colour content of the qqg system 191]. 
The einbein fields are to be found from the constraint conditions |l6| : 

dUg dflg dflg 

Thus, to quantise the system, one should find the einbeins from Eq. ([71) and substitute them 
back to the Hamiltonian (I5|). In such a way, einbeins would become entangled functions 
of coordinates and momenta. To avoid these complications, an approximate einbein field 
method is used in the QCD string model calculations: einbeins are treated as c-number vari- 
ational parameters. The eigenvalues of the spinless Hamiltonian ([5]) are found as functions 
of fig, fig, and fig and, finally, minimised with respect to the einbeins. With such simpli- 
fying assumptions the spinless Hamiltonian Hq takes an apparently nonrelativistic form, 
with einbein fields playing the role of the constituent masses of the quarks and the gluon. 
These quantities, however, are not introduced as model parameters, but are calculated in a 
relativistic formalism. Indeed, the procedure of taking extrema in the einbeins is nothing 
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but the summation of the entire series of relativistic corrections to the leading would-be 
nonrelativistic eigenenergies. 

The einbein field method allows one to estimate corrections to the leading potential regime 
([5]). First, as was discussed before, one takes into account the contribution of the string 
inertia, expanding the square roots in the expression ([3]) to the first order in 2 — this gives 
the string correction Kjtr- Second, employing the Feynman-Schwinger representation for the 
Green's functions of spinning quarks and gluons one can extract the spin-dependent part of 
interaction, as described in Ref. [20|, hereafter denoted as Vsd- Finally, the FCM method 



allows one to calculate the nonperturbative selfenergy of the quarks which, as was shown 



in Ref. 2l|, provides an overall shift of the hadron mass, as required by phenomenology. 
We use the notation C for this contribution. Thus the full form of the hybrid Hamiltonian 
reads: 

H = Ho + V^tr + + C. (8) 

To specify extra terms in Eq. ([8]) let us first introduce various angular momentum oper- 
ators: 

Llg = [fsiPg], L2g = - [r23Pg], ^Ig = -[r^lPg], ^2? = [^23^9], Lg = [fuPg], Lg = -[f^Pg], 

(9) 

where 

ri2 = rg- rg, fgi = fg - fg, ^23 = Tg - Vg. (10) 

In terms of these coordinates and angular momenta the string correction Kjtr is trivially 
calculated from Eq. (IHl) and takes the form (due to the symmetry of the problem we set 

6ri3 I ^2 ^2 I yq^gl- 



The derivation of the spin-dependent potential can be found in Ref. 20|], with the result: 

^SD = Vif + Vii> + Vss + V^f + V^^\ (12) 

where the superscript LS stands for the spin-orbit interaction, SS — for the hyperfine 
interaction, ST — for the spin-tensor, and 



^Ls'^ = -ji^ ( - L,sg ) , (13) 



4/i2r 
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2ri3 1^ /i2 ^2 j 2rl^\2^ig fij fig 2rf^ \2fi fig) /i \q^q)' 

(14) 

Anas ^ ^ c-/^ X 47ra5 ^ ^ j.,^ , Anag ^ ^ \ ^\ 
»/ss = — -^SqSqd[r 12) H SqSgd[ri^) H SgSgd^ras), (15) 

^st'^ = -iSt^ (3(3/12) (Sgri2) - r?25;s5) , (16) 

^s? = mn,)isgn,) - rlsqSg) + (J ^ . (17) 

Formally, these spin-dependent terms coincide with the well-known Eichten-Feinberg- 
Gromes ones 22j. Notice, however, that the Eichten-Feinberg-Gromes spin- dependent po- 
tential comes out from an expansion of the interaction in the inverse powers of quark masses, 
whereas the potential Vsd derived above contains effective constituent masses fi and fig in the 
denominators instead of current masses. Due to confinement, these constituent-like masses 
are always large, of order of the confinement scale (~ ^/o' ~ 400 MeV) or larger, even for 
massless particles. So the result for Vsd is applicable to the case of light quark flavours, as 
well as to the case of massless gluons. The only approximation made in order to derive the 
spin-dependent potentials ( fT3i) -( |T7l) is the Gaussian approximation for field correlators — 
see Ref. j^] for the discussion. 

Finally, the constant C is 

C = -—rf{m,Tg), (18) 
fin 

where ri{m,Tg) is a universal function of the quark mass m and the gluonic correlation 



length Tg — see Ref^ 
formalism fin Ref. 



21 



2l| for the explicit form of this function and for further details of the 
] the gluonic correlation length is denoted as 6). 



III. THE SPECTRUM OF HYBRIDS 

The form ([5]) allows one to separate the centre-of-mass motion in a standard way, intro- 
ducing the Jacobi coordinates as 

r = rg- Tg, p = rg -— - — = rg -— . (19) 

M'q 1 flq ^ 

In terms of these Jacobi variables the Hamiltonian Hq in the centre-of-mass frame can 
be written as 

H^ = — +fi+^ + ^ + -9— + ^rsi + ar23 + VcouU (20) 

fl 2 2/ii2 2/ii2,3 
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where 



/ii2 = f, ^3 = ^, M = 2^^ + ^^,, (21) 



and p and Q are the momenta conjugated to the Jacobi coordinates r and p, respectively. 

First, we note that the relative angular momenta Iqg and Ig are not conserved in the 
three-body qqg system, though the requirement of a given 

jPC 

imposes restrictions on their 

possible values. On the other hand, the zero-order Hamiltonian (1201) conserves the total 
quark spin Sgg. In accordance with Eqs. ([1]) and ([2]), states with magnetic and electric 
gluons have different total quark spin, so they are not mixed in the leading order. Thus one 
can employ the well-known hyperspherical formalism to calculate the zero-order spectrum 
and w.f.. 

As we are interested in hybrids with a magnetic gluon, we use trial w.f.'s of the form 

= Hr,p)Soiqq)Y,ClZu2pYiuAp)Si.M, (22) 

for the vector hybrid, and 

I J-+)„ = $(r,p) C(;:,,,S^,,iqq)Y,Cl!::i.,pY,.riP^ (23) 

for its siblings. Here Si^{g) is the spin w.f. of the gluon, So{qq) and Siy{qq) are the singlet 
and triplet spin w.f.'s of the qq pair. The "radial" w.f. $(r, p) depends on its arguments in 
the form of the hyperspherical radius i?, 

^ + ^-f^^. (24) 

Necessary formulae of hyperspherical formalism in three-body systems can be found in 
Appendix |A] (see also Ref. j24| for more details). In actual calculations the radial w.f. was 
chosen in the Gaussian form, 

$(r,p)=exp(^-^/?2Mi?2^, (25) 
with (5 being the variational parameter. Then the eigenvalue Mq is given by 

Mo = h{po,pgo,Po), (26) 

where 

M/i,/i„/5) = ^^^p, ^{f,p)=pYUmr,p). (27) 
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Parameter m, GeV a, GeV^ ag rj 

Value 1.48 0.16 0.55 0.29 



TABLE I: The set of parameters used for the numerical evaluation. 



State 


^^0 




'Pi 


'Pi 


3p 

Pq 


'P2 


Exp 


2.980 


3.097 


3.526 


3.511 


3.415 


3.556 


Theor 


2.981 


3.104 


3.528 


3.514 


3.449 


3.552 



TABLE IL Masses of S- and P-level charmonia, in GeV. 

and /io, /igO; and (3q come as a selfconsistent solution of the set of three coupled equations, 

dh ^ dh ^ dh ^ 
df5 ' dfi ' dfig 

The model parameters {m, a, Og, r]} are fixed by evaluating the spectrum of conventional 
cc IS"- and IP- states within the same formalism (in Appendix [B] we give necessary details 
of calculations for the conventional charmonium). The set of parameters is given in Tabled 
and the results for charmonia levels are listed in Table HTl The value of ?7 = 0.29 corresponds 
to the charmed quark mass given in Tableland to the gluonic correlation length ~ 0.2 
fm, which complies well with the one measured on the lattice (around 0.2 ^ 0.3 fm M). 

With these parameters the variational procedure described above yields 

Mo = 4.573 GeV, (29) 

for the zero-order hybrid mass, while the extremal values of the effective masses are 

/io = 1.598 GeV, fXgo = 1.085 GeV. (30) 

Notice that the gluon effective mass figo appears to be rather large, and of the same order 
of magnitude as the effective charmed quark mass fiQ. 

All corrections to the leading regime (l29l) are calculated as perturbations, with the sub- 
stitution fi ^ fiQ, fig ^ figQ, (3 I3q. There are two types of such corrections: one that does 
not depend on quark spin and the other, which depends. The former correction provides an 
overall shift with respect to zero-order regime (l29l) . The spin-dependent correction removes 
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the degeneracy between the four states (122!) and (123!) and, in principle, is responsible for the 
mixing of the magnetic gluon states with electric gluon ones. In what follows we neglect 
such a mixing. 

The simplest correction of the first type is the selfenergy correction which shifts the 
zero-order hybrid mass downwards, 

AM,elfenergy = "28 MeV. (31) 

The string correction does not depend on spins either, and it is calculated to be 

AMstring = -52 MeV. (32) 

There is also a mass shift due to the gluon spin-orbit force, common for all four states 
( I22!) and (l23l) . which comes from the terms in (lT4l) proportional to the operator Sg. This 
yields: 

AM^g = -103 MeV. (33) 

Finally, to calculate the spin splittings, it is convenient to rewrite the operators Sq and 
Sq in terms of operators Sqq and S, 

Sqq=\{.Sq + Sq), t = ]^{Sq-Sq). (34) 

We notice then that, once the operator S is antisymmetric with respect to the permutation 
g ^ g, it flips the spin of the qq pair and, as such, is neglected in our calculations. Then, 
after tedious but straightforward calculations, one arrives at the following expressions for 
spin splittings: 



AM(r-) = 


, 3 A (gq) 




AM(0"+) = 




- 2A, 


AM(1-+) = 


1a (99) 


-A, 


AM(2-+) = 


1a (99) 

4^ss 


+ A, 



(35) 



where 

A^f ) = 9 MeV (36) 
comes from the quark-antiquark spin-spin interaction, and 

A = A(f + Ag + AsT. (37) 
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0-+ 


1-+ 


1"- 


2-+ 


Mass 


4.252 


4.320 


4.397 


4.457 



TABLE III: Masses of charmonium hybrids, in GeV. 

Individual contributions in Eq. fl37|) are 

a1?|^ = 9 MeV, = 24 MeV, Ast = 35 MeV, (38) 

coming from the quark-gluon spin-spin interaction, the spin-orbit interaction proportional 
to the quark and the antiquark spin, and the spin-tensor interaction, respectively. These all 
together give 

A = 68 MeV. (39) 

The ultimate numerical results for the hybrid masses are given in Table UTTt spin splittings 
are established to yield: 

M(0"+) < M(l"+) < M{r~) < M(2-+). (40) 



IV. COMPARISON TO OTHER APPROACHES AND LATTICE 

CALCULATIONS 



The story of hybrid meson studies started with the bag model calculations 26j, where 
the lowest charmonium hybrid mass of about 4 GeV was obtained, and the splitting pattern 
(140|) was found. 

In the flux tube model the hybrid excitations are visualised as phonon-type exci- 
tations of the string connecting the quark-antiquark pair, and a certain correspondence is 



established in Ref. [28j between the excited flux tube and the constituent gluon approaches. 
Eight lowest ccg hybrids are predicted [29] to reside around 4.1-f-4.2 GeV, and 1 and J ^ 
states are among those. However, there is a discrepancy in quantum numbers of the flux 
tube hybrids and the ones with a constituent gluon: the constituent gluon carries colour 
and spin. As a result, the P-even flux tube hybrids have 

j^t: = 0+-,l+-,2+-,l++ (41) 
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quantum numbers, while in constituent gluon models these are 



(42) 



hybrids with the electric gluon. Spin splittings of the flux tube hybrids due to the long-range 



Thomas precession were calculated in Ref. 



30| to be small, while the splittings reported in 



the present paper are much larger, and come mostly from perturbative short-ranged forces. 

The constituent gluon model with pairwise forces was presented in the pioneering work 
in Ref. j^, and in Ref. jsil. The mixing between magnetic and electric gluon hybrids 
was calculated in Ref. {sil, with rather controversial results: in the first entry of Ref. 31 1 
the mixing was found to be small, while in the second entry it is claimed to be substantial 
(though the details are not given there). The unmixed states are found in the same mass 
region as in the present work. 

There exist results 32] of the QCD string model calculations in the einbein field for- 
malism. The l""*" charmonium hybrid mass was found to be 4.2 ± 0.2 GeV. However, this 
result cannot be directly compared with ours, as the adiabatic approximation for quarks is 
employed there. 

Finally, in a potential model, with the cc pair considered as a colour-octet source, the 
tensor hybrid was predicted at 4.12 GeV [3^. 

Lattice simulations deal mostly with exotic quantum numbers, with l""*" charmonium 



hybrid residing at about 4.4 GeV — see Ref. 
recent results we would like to mention Ref. 
charmonium hybrid, and Ref. 



34j | and references therein. Among more 
3g, which gives 4.405 ±0.038 GeV for the 1"+ 



36j where, for the first time, a 1 state was found, excited 



by the hybrid meson operator, with the mass of 4.379 ± 0.149 GeV (the authors claim that 
they have found a radially excited vector hybrid; such an interpretation was criticised in 
Ref. 33] ). As seen from our Table UTTl the agreement of lattice results and our findings is 
quite good. 



V. HYBRID DECAYS 

The selection rule forbidding the decay of 1 and J ^ hybrids with magnetic gluon into 
the final states was established for the constituent glue model in Refs. 0, Q]. 

As the decay takes place via gluon quark- antiquark pair transition, the amplitude is 
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DDo 


^♦Do DDA^Pi) 


D*Di(^Pi) 


DDii^Pi) 


D*Di(^Pi) 


DD2 


D*D2 


1— 


1 


1 
1 


1 
1 


1 
1 




1 /5 


V ^ 


~2 


2 


2v^ 




-2V6 






1 












72 










1-+ 


1 1 


1 


1 


1 








2^/2 




4 






2-+ 




1 




3 


1 






2^/2 




4 


2^2 


4 



TABLE IV: Spin-recoupling coefficients for the hybrid states listed in Table IIIII Here D^*') is an 
S'-wave D''*'*-meson and Dj is a P-wave Z)-meson with the total momentum J. A proper charge 
conjugation is implied. 

proportional to the overlap of the initial and final state w.f.'s with the pair creation operator. 
Thus it vanishes if the final state mesons have identical spatial w.f.'s. The latter is definitely 
a eood approximation in the case of D^*^D^*\ Similar symmetry considerations hold true in 
the t„be .odel and .ee,n to be .atKe. ,ene.a, a, demonstrated in Ref. «. 

ection rule based on the spin content of quarks in the final state, 



There exists another se^ 
established in Refs. jsl, IsS] and, quite recently, in Ref. [39j. Assuming that i) the spin of 
the initial qq pair does not flip in the decay, and ii) the qq pair created in the decay is in 
spin-triplet, one can define the relative strength of matrix elements for decays into various 
final meson pairs. This selection rule is rather powerful: for example, for a vector hybrid, 
the initial quark pair is in spin-singlet, while for a conventional vector quarkonium it is 
in spin-triplet, with clear discrimination between two possibilities. Thus, measuring the 
relative rates of various S- and P meson pairs, one could distinguish between a vector 



quarkonium and a vector hybrid 39|, |40| . 



For hybrid decays under consideration, only S-wave amplitudes are of immediate rele- 
vance (the D-wave ones are suppressed because of a limited phase space), and the corre- 
sponding spin-recoupling coefficients are given in Table HV] (decay rates are proportional to 
the squares of these coefficients). The coefficients exhibit a sum rule: if the masses and w.f.'s 
of the initial and final states involved could be taken identical, the total effect of coupling 
to open charm mesons is identical within the whole hybrid multiplet (though individual 
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contributions from various D-meson channels differ). A similar sum rule holds true for the 



conventional cc charmonium couplings to D-mesons, as found in Refs. 



41 



42] and formu- 



lated in Ref. 43| as a general theorem. Most straightforward consequence of this sum rule 
is the following observation: the effects of mesonic loops over the spectra could be quite 
large, but numerically they are similar for all low-lying charmonia. For hybrids, however, 
the latter is not the case as, in accordance with the results of Table llllt different members 
of hybrid multiplet reside among different thresholds. 

VI. DISCUSSION 

First, we notice that both lattice calculations and our findings place the vector hybrid at 
4.4 GeV, substantially higher than the y(4260). One should have in mind, however, that 
the lattice result 36| for a vector hybrid comes with a large error, and the accuracy of the 
einbein field method is not better than 5% in the binding energy, as estimated in Ref. 11| . 
so the discrepancy could appear to be not very significant. On the other hand, the influence 
of open charm channels is not taken into account in the present approach, and the said 
influence could be large. 

Indeed, as seen from Table IIVI there is a significant coupling of the vector hybrid to the 
DDii^Pi) channel. There are two Di-mesons, a narrow one with the mass of 2420 MeV and 
the width of 20 MeV, and a broad one with the mass of 2430 MeV and the width of 380 
MeV 4J], which are (unknown) mixtures of the ^Pi and ^Pi states. Thus the vector hybrid 
state should be attracted to the corresponding thresholds, which are tantalizingly close to 
the measured mass of the y(4260). 

It is interesting to mention in this regard another enigmatic vector state (or, maybe, even 



two states!), namely y(4325) from BaBar 45| and y(4360) from Belle 46|], both seen in 



the initial state radiation process, with the masses consistent with each other and with the 
width of the Belle state being two times smaller than that of the BaBar one. The masses of 
these new y's are close to another relevant threshold, the D*Dq one (we assume, following 
the Belle paper [4:7], the scalar D-meson to be at 2308 MeV, with the width of about 270 
MeV). So the coupling of the hybrid vector to DDi and D*Dq could be responsible for the 
formation of two near-threshold states. 

The case of the pseudoscalar hybrid is simpler: in accordance with Table IIVI half of the 
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decay strength goes to the DDq channel with the nominal threshold at 4.18 GeV, so this 
hybrid would feed the structure in the DDtt final state, with this structure being broad due 
to large Dq width. 

Exotic hybrid is estimated to be 80 MeV lighter than the vector one, with the mass very 
close to the D*Dq threshold. As seen from Table IIVI the coupling to D*Dq is two times 
larger than the one for the vector case, so, as the i5o(2308) is very broad, this hybrid has 
more chances to disappear in the D*D7t continuum. 

As to the tensor hybrid, with the bare mass of 4.457 GeV, it could survive as a resonance, 
because the most prominent channel, D*Di, opens nominally only at 4.43 GeV. 

VII. CONCLUSIONS 

In this paper we calculated the masses of low-lying charmonium hybrids with magnetic 
gluon in the framework of the Field Correlator Method for QCD. The QCD string approach 
is employed to estimate spin-dependent corrections for the 1 , ^, 1^"'', and 2 ^ hybrid 
states. The spectrum is calculated without fitting parameters, as all the parameters of the 
effective Hamiltonian are fixed by reproducing cc charmonium levels. Our results are in 
good agreement with lattice data. 

Decay modes of hybrids are investigated and predictions are made for the relative rates 
of the decays into S- and P-wave D-mesons for all four states of the lowest hybrid multiplet. 

The calculated mass of the vector hybrid is 4.397 GeV, substantially higher than the mass 
of a promising hybrid candidate y(4260). We argue that strong coupling of the vector hybrid 
to the DDi and D*Dq modes can cause considerable threshold attraction, making vector 
hybrid bare state responsible for the formation of near-threshold y(4260) and y(4325) 
states. 
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APPENDIX A: SOME DETAILS OF THREE BODY KINEMATICS 

In a three-body system of particles with masses /ii, ^2, /^s three sets of Jacobi coordinates 
can be defined. One is 

ri2 = n - ^2, P3 = r3 ■ , (A.l) 

/ii +1^2 

and others {f3i,Q2 and f23,Qi) are obtained from it by cyclic permutations of the particle 
indices. The hyperspherical radius 

where 

yUl + /i2 M 

is invariant under such cyclic permutations. Similarly, 

d^n2d^P3 = d^r3id^p2 = c/VasC^Vi- (A-4) 
Angular momenta can be conveniently represented as 

Lig = fsi X ^i^rgi X Q2, Lig = fgi X H ^i^fai x Q2, 

P + Pg P + Pg 

L2g = r23 X P23 H ^^^r23 X Qi, L2g = r23 X P23 ^;^r23 X Qi, (A. 5) 

1 - - 1 - 

Lq = fi2 X pi2 - -fi2 X (^3, Lq = fu X pi2 + -fu X Q3, 

where the momenta pij and Qa; conjugated to corresponding Jacobi coordinates are intro- 
duced, and pi = P2 = P, P3 = Pg is substituted. 
The formula 

is extensively used in calculations of various matrix elements. 
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APPENDIX B: SPECTRUM OF CONVENTIONAL CHARMONIA 

In this Appendix we give some details of evaluation of the spectrum of conventional 
charmonia which was used in order to fix the set of model parameters. Further details of 



calculations for conventional mesons can be found in Refs. 



11 



12|. 



The spinless centre-of-mass Hamiltonian of the charmonium reads: 



i/f) = 2Vp2 + m2 + ar-^^, (B.l) 
where r is the interquark distance. The einbein fi is introduced then, 

2^/7^^^!!^ + ^ + ^. (B.2) 

The variational procedure for the charmonium is similar to the one described in the text 
body for the hybrid, that is, we take a Gaussian trial w.f., 

^('=2)(r) = exp (^—fip^r'^^ , (B.3) 

and evaluate the average 

h^'^Xfi^P) = ^ . (B.4) 

yi^^ / /\[/(cc) Kl/(cc)\ ^ ^ 

Then, taking extrema in the variational parameters fi and (3, we find the "bare" charmonium 
mass and the actual value /io to be substituted to the corrections to the Hamiltonian fIB.ip . 



As in case of the hybrid, these corrections can be classified as spin-independent and spin- 
dependent. 



The spin-independent corrections to the Hamiltonian (IB.ip are the string correction 18| , 



= (B.5) 

6 fir 

with L being the quark-ant iquark angular momentum, and the selfenergy correction C which 
coincides with the one for the hybrid, given by Eq. ( |T8l) (this is due to the fact that the 



quark content is the same for both the charmonium and charmonium hybrid 2l|). 
The spin-dependent corrections are: 

vS' = vif + vff' + yP^ (B.6) 

where 

V^f = -:^LS + ^^LS, (B.7) 
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a(cc) 
SS 



327ias 



(B.8) 



9/^2 




(B.9) 



In Table HTl we compare the predicted masses of the IS"- and IP-wave charmonium states 
with the experimental data. The set of parameters used in the calculations is given in 
Table [H We find a good agreement of our predictions with the data and use the same set of 
parameters to evaluate the masses of the hybrids. 
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